Tensorial slip of super-hydrophobic channels 
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We describe a generalization of the tensorial slip boundary condition, originally justified for a thick 
(compared to texture period) channel, to any channel thickness. The eigenvalues of the effective 
slip length tensor, however, in general case become dependent on the gap and cannot be viewed as 
a local property of the surface, being a global characteristic of the channel. To illustrate the use 
of the tensor formalism we develop a semi-analytical theory of an effective slip in a parallel-plate 
channel with one super-hydrophobic striped and one hydrophilic surface. Our approach is valid for 
any local slip at the gas sectors and an arbitrary distance between the plates, ranging from a thick 
to a thin channel. We then present results of lattice Boltzmann simulations to validate the analysis. 
Our results may be useful for extracting effective slip tensors from global measurements, such as 
the permeability of a channel, in experiments or simulations. 

PACS numbers: 47.ll.-j, 83.50.Rp, 47.61.-k 



I. INTRODUCTION 

With recent advances in microfluidics [ll, 01 , renewed 
interest has emerged in quantifying the effects of sur- 
face heterogeneities with different local hydrophobicity 
(characterized by a local scalar slip [3,3), on fluid mo- 
tion. In this situation it is advantageous to construct 
the effective slip boundary condition, which is applied at 
the hypothetical smooth homogeneously slipping surface, 
and mimics the actual one along the true heterogeneous 
surface 0, Q . Such an effective condition fully character- 
izes the flow at the real surface and can be used to solve 
complex hydrodynamic problems without tedious calcu- 
lations. Well-known examples of such a heterogeneous 
system include super-hydrophobic Cassie (SH) surfaces, 
where trapped gas is stabilized with a rough wall tex- 
ture, leading to a number of 'super' properties, such as 
extreme non- wettability and low hysteresis ■ For these 
surfaces effective slip lengths are often very l arg e (8l-[ll| 
compared a smooth hydrophobic coating p^ - ll7| . which 
can greatly reduce the viscous drag and impact transport 
phenomena in microchannels [H . 

The concept of effective slip was mostly exploited for 
thick (coinpared to the texture characteristic length, L) 
channels [l|, [l^ , where for an anisotropic texture it was 
shown to depend on the direction of the flow and is a 
tensor [l^, boff = represented by a symmetric, 

positive definite 2x2 matrix 



diagonalizcd by a rotation with angle Q 

„ / COS0 sinG 
^® " I -sine cose 
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For all anisotropic surfaces its eigenvalues and 6^ 
correspond to the fastest (greatest forward slip) a.nd slow- 
est (least forward slip) orthogonal directions [l^. In the 
general case of any direction e, this means that the flow 
past such surfaces becomes misaligned with the driving 
force. This tensorial slip approach, based on a consid- 
eration of a 'macroscale' fluid motion instead of solv- 
ing hydrodynamic equations at the scale of the individ- 
ual pattern, was supported by statistical diffusion ar- 
guments [13, and was recently justiflcd for the case of 
Stokes flow over a broad class of periodic surfaces Q. 
Note that an effective slip in a thick channel situation is 
a characteristics of a heterogeneous interface solely (be- 
ing expressed through its parameters, such as local slip 
lengths, fractions of phases, and a texture period) [5l.l20l|. 

It was however recently recognized and justified by us- 
ing the theory of heterogeneous porous materials [2l| that 
a similar concept of effective slip can be also exploited for 
a flow conducted in a thin channel with two conflning sur- 
faces separated by a distance H L. In such a situation 
a natural definition of the effective slip length could be 
based on a permeability of a hypothetical uniform chan- 
nel with the same flow rate. An effective tensorial slip 
is then determined by flow at the scale of the channel 
width, and depends on H [2l|. This points to the fact 
that an effective boundary condition reflects not just pa- 
rameters of the liquid-solid interface, but also depends 
on the flow configuration 

The power of the effective slip approach and the super- 
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FIG. 1: Sketch of the SH stripes (a): B = 7r/2 corresponds 
to transverse stripes, whereas O = to longitudinal stripes; 
(b) situation in (a) is approximated by a periodic cell of size 
L, with equivalent flow boundary conditions on the gas-liquid 
and solid-liquid interfaces. 



II. MODEL AND GOVERNING EQUATIONS 

The basic assumptions of our theoretical model are as 
follows. We consider a channel consisting of two paral- 
lel walls located at y = and y = H and unbounded 
in the x and z directions as sketched in Fig[TJ The up- 
per plate represents a no-slip hydrophilic surface, and the 
lower plate is a SH surface. The origin of coordinates is 
placed it the plane of a liquid-gas interface at the center 
of the gas sector. The x axis is defined along the pressure 
gradient. This (SH vs. hydrophilic) geometry of configu- 
ration is relevant for various setups, where the alignment 
of opposite textures is inconvenient or difficult. Besides 
that the advantage of such a geometry is that it allows to 
avoid the gas bridging and long-range attractive capillary 
forces [25| , which appear when we deal with interactions 
of two hydrophobic solids [2^, [27l . 

As in previous publications j20l [23l . [28j , we model the 
SH plate as a flat interface with no meniscus curvature, 
so that the SH surface appears as perfectly smooth with 
a pattern of boundary conditions. The latter are taken 
as no-slip (bi = 0) over solid/liquid areas and as par- 
tial slip (&2 — b) over gas/liquid regions. We denote as 
S a typical length scale of gas/liquid areas. The frac- 
tion of the solid/liquid areas is denoted = (L — S)/L, 
and of the gas/liquid areas ^2 = 1 — = (^/i. In this 
idealization, by assuming a flat interface, we neglect an 
additional mechanism for a dissipation connected with 
the meniscus curvature [H, [2^ Ho] . 

The flow is governed by the Stokes equations 



lubrication potential of SH surfaces were already illus- 
trated by discussing several applications. In particular, it 
has been shown that optimized SH textures niay be suc- 
cessfully used in passive microfluidic mixing 0, 22 1 and 
for a reduction of a hydrodynamic drag force 23 , [l^l , 
and that the effective slip formalism represents a useful 
tool to quantify properties of SH surfaces in thick and 
thin channels. In many situations however the dramatic 
changes in flow happen when the two length scales H and 
L arc of the same order. In this work we generalize the 
deflnition of the effective slip length tensor Eq. ([1]) to an 
arbitrary channel thickness and validate this approach by 
means of lattice Boltzmann (LB) simulations. 



The structure of this manuscript is as follows. In sec- 
tion II we describe our model and formulate governing 
equations. In section HI we analyze a tensorial perme- 
ability of the parallel-plate channel with one (anisotropic) 
SH and one hydrophilic surface and give some general 
arguments showing that the tensorial relation, Eq. ([1}, 
should be valid at arbitrary H/L. Section IV contains 
semi-analytical results for a striped SH surface, and in 
section V we describe our LB simulation approach. In 
section VI simulation results are presented to validate 
the predictions of the tensorial theory and to test our 
analytical results for the asymptotic limits of thick and 
thin channels. 



TyV^u = Vp, 



V-u = 0, 



(3) 



(4) 



where u is the velocity vector, and the average pressure 
gradient is always aligned with the x-axis direction: 



(Vp) = (-a,0,0) 



(5) 



The local slip boundary conditions at the walls are de- 
fined as 

u{x,0, z) ^ b{x,z) ■ —{x,0,z), y • u(a;,0,z) = 0, (6) 

oy 



u{x,H,z) ~ 0, y ■ u{x,H,z) = 0. 



(7) 



Here the local slip length b{x, z) at the SH surface is 
generally the function of both lateral coordinates. 

We intend to evaluate the effective slip length 6off at 
the SH surface, which is as usual defined as 



&off 



my 

where (. . .) means the average value in the plane xOz 



(8) 
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III. GENERAL CONSIDERATION 

In this section, we evaluate the pressure-driven flow in 
the parallel-plate channel with one anisotropic SH and 
one hydrophilic surface and investigate the consequences 
of anisotropy. For an anisotropic texture, there are sev- 
eral possible ways to define effective slip lengths. A nat- 
ural approach is to define a slip-length tensor via 




(9) 



by analogy with a hypothetical uniform channel. Such 
a definition was earlier justified for thin SH channels, by 
using the lubrication limit [2l|. Below wc argue that 
the same result is obtained for a channel of an arbitrary 
thickness with a tensorial generalization of the Navier 
boundary condition, where beff is a global measure of 
the effective slippage of the channel. 

For a mathematical justification of the above state- 
ment we first rewrite Eqs. ([3]) and ([5]) as 



TjV^U) = (Vp), {Vp) = -aSi, 



(10) 



where (U) is the depth-averaged velocity. The two 
boundary conditions that apply at the channel walls can 
then be written as follows: (U) = at the upper surface, 
and uniform tensorial slip, 



(u) = beff -(ayU), 



(11) 



at the lower surface, with beff defined according to 
Eqs. (dl-©. The solution for the depth-averaged ve- 
locity is then given by Q 



with 



2i] 2i] 



2,; 



Hbl^ cos2 e + Hb^ff sin^ e + 



(13) 



cff^eff 



iH + bl^){H + b^^) 

H{blg - fe^g) sine cose 
{H + bl)iH + b^^) 



In linear response, the averaged fiow rate, (Q), is pro- 
portional to (Vp) via the permeability tensor, k: 



(Q) 



-k-(Vp). 



(14) 



Integrating the velocity profile across the channel we ob- 
tain 



H 



(Q) 



(U(2/)) dy, 



(15) 



with the components 

j3 



?? 12 



77 4 



(16) 



The latter may be rewritten as 



(Q)^ = - (^/cllcos^e + fc-Lsin^ej , (17) 

(g)^ = ^ (^/c" - /c-L^ sine cose, (18) 

provided the two tensors, k and bed, are coaxial and the 
rigorous relationship between their eigenvalues is given 
by Eq. ©. This unambiguously indicates that the two 
definitions of the slip length are equivalent. It becomes 
also apparent that Eq. ^ implies that generally 
depends on the separation H . 

Note some similarity to a prior work [2l|, [3l|. The 
current consideration, however, is valid for arbitrary H / L 
and b/L, including the limit of thick channels, where boff 
becomes a local property of the surface on scales much 
larger than the texture characteristic length. 

Consider now a situation where the 'fast' axis of great- 
est forward slip of anisotropic texture is inclined at an 
angle e to the pressure gradient. This problem can be 
solved explicitly as follows. The downstream effective 
permeability of the channel can be expressed in terms of 
the effective downstream slip length as 



(19) 



Following [19[, it can also be obtained from the perme- 
ability tensor: 




fcllfcj 



fell sin^ e + fc-L cos^ e 



(20) 



By substituting Eq. ^ into Eq. ((20| and after subtract- 
ing the latter from Eq. (|19p we express an effective down- 
stream slip length in the form 



b^^H + Abl^bi^ 



\^-bi^)H cos' Q ^21) 



H + Ab\-A{bl^-b^^)cos'Q 



Note that in the general case b),^ depends on H and 

&||^(-ff). For this reason, 6^^'' cannot be viewed as a local 
property of the SH surface, except as in the thick channel 
limit. Instead, it is generally the effective slip length of 
the SH channel and thus its global characteristic. 

Finally, we emphasize the generality of Eqs. (f20| and 
(|2ip , which follow only from the symmetry of the effective 
permeability and slip- length tensors for linear response. 
Similar formulae have been obtained before in a few par- 
ticular calculations [HI |32|, but the present derivation 
is valid regardless of the thickness of the channel and is 
independent of the details of the textured surface. There 
could be arbitrary patterns of local slip lengths, and the 
latter could itself be a spatially varying tensor, reficcting 
surface anisotropy at a smaller (possibly atomic) scales. 
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IV. THEORY FOR STRIPED PATTERNS 

To illustrate the general theory, in this section we 
focus on flat patterned SH surfaces consisting of peri- 
odic stripes, where the local (scalar) slip length b varies 
only in one direction. The problem of flow past striped 
SH surfaces has previously been studied in the con- 
text of a reduction of pressure-driven forward flow in 
thick [2^, m, H^l and thin f2li channels, and it is di- 
rectly relevant for mixing [sl |22| | , a. nd a generation of a 
tensorial electro-osmotic flow [5|, [Sa, |3^ . Here we elabo- 
rate on a previously published anzatz [5| and present the 
theory for an arbitrary gap, which in the asymptotic lim- 
its describes situations of thin and thick channels. The 
mathematical anal ysis we use here is similar to a tech- 
nique exploited in |34l | for a thick channel configuration. 
The crucial difference with Ref. [s^], however, is that we 
consider an arbitrary gap, which means that the effective 
slip is a function of the channel thickness as discussed 
above. 

For transverse stripes, we have u = {u{x, y),v{x, y), 0), 
w(x,0) = b{x)uy{x,0), and v{x,0) = 0. For longitu- 
dinal stripes, the flow is also two dimensional: u — 
{u{y, z),v{y, z),0), u{0,z) = b{z)uy{0, z), and v{0,z) — 
0. As the problem is linear in u, we seek a solution of 
the form 



The Fourier method yields a general solution to this prob- 
lem: 

u,iy,z) = (Moy + Po) (27) 

00 

+ ^(A/„e^"^ + Pne-^"«) cos(A„z), 

with A„ = 2'Kn/L. The sine terms vanish due to symme- 
try. Condition ([25)) leads to 



ui{y,z) 



OG 

^P„cos(A„z)e-^"''(l 



(28) 



n=l 



Applying the boundary condition (|26p we then obtain a 
trigonometric dual series: 



1 + — ) + ttn [1 + bXn COth(A„g)] COs(AnZ) 



b—, < z< (5/2, 
21] 



(29) 



ao 



^ a„ cos(A„z) = 0, S/2<z<L/2, (30) 



U = Uo + Ui, 



(22) 



where Uq is the velocity of the usual no-slip parabolic 
Poiscuille flow 

uo = (wo,0,0), uo = -^y^ + ^y (23) 

2i] 2i] 

and Ui is the SH slip-driven superimposed flow. 



A. Longitudinal stripes 

In this situation the problem is homogeneous in x- 
direction {d/dx = 0). The slip length b{x,z) — b{z) 
is periodic in z with period L. The elementary cell is 
determined as b{z) = 6 at \z\ < 6/2, and b{z) = at 
S/2 < \z\ < L/2. In this case the velocity Ui ~ (mi,0, 0) 
has only one nonzero component, which can be deter- 
mined by solving the Laplace equation 



V^ui{y,z)^0, 



(24) 



with the following boundary conditions deflned in the 
usual way as 



wi(i7,z) =0, 



ui(0,z) = fo(z) 



/ aH 


dui 






dy 





(25) 



(26) 



where 



-2X„H 



n>l. 



Dual series (|29l). (|30l) provide a complete description of 
hydrodynamic flow and effective slip in the longitudinal 
direction, given all the stated assumptions. These equa- 
tions can be solved numerically (sec Appendix E| . but 
exact results are possible in the limits of thin and thick 
channels. 

For a thin channel, H <C L, we can use that 
cothf|j_^g = + 0{t). By substituting this expres- 
sion into and and keeping only values of the first 
non- vanishing order [ssj . we find 



s/2 

2 r aH 



ao 



whence |2lf 



.dz^^Ji!^ (31) 



2i] 1 + b/H 2ri H + b 



"eft 



H 



(32) 



This is an exact solution, representing a rigorous upper 
Wiener bound on the effective slip over all possible two- 
phase patterns in a thin channel. In order to gain a sim- 
ple physical understanding of this result and to facilitate 
the analysis below, it is instructive to mention the two 
limits d ^ that follow from Eq. When H <^b,L 

we deduce 



~ —H oc H, 



(33) 
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and when 6 ^ iJ <C i we get a surface averaged slip 



^off 



oc b. 



(34) 



In the limit of a thick channel, 7? 3> we can use 
that coth(t — !■ oo) 1 and the dual series ([^ -([50 )) can 
be solved exactly to obtainQ 



^cff 



L 

TT 



In 


sec 










L 

1 + — In 

TTO 


sec 


it) 


+ tan 


it) 





(35) 



This expression for an effective slip length depends 
strongly on a texture period L. When 6 <C L we again 
derive the area-averaged slip length 



^off 



oc b. 



When b ^ L, expression takes the form 
L 



(36) 



In 



L^bM TT 



2 



that coincides with an earlier result 
perfect slip (6 oo) case. 



oc L, (37) 
obtained for a 



B. Transverse stripes 

In this case it is convenient to introduce a stream 
function tp{x,y) and the vorticity vector u}{x,y). The 
two-dimensional velocity field corresponding to the 
transverse configuration is represented by u(x, y) = 
{dtjj/dy,—dtp/dx,0), and the vorticity vector, Lj{x,y) = 
Vxu = (0,0, cj), has only one nonzero component, which 
equals to 



OJ 



(38) 



The solution can then be presented as the sum of the base 
fiow with homogeneous no-slip condition and its pertur- 
bation caused by the presence of stripes as 



if.' = 'i'o + V'l : uj ^ flo + uji, 



(39) 



where '^q and ^lo correspond to the typical Poiscuille flow 
in a flat channel with no-slip walls: 



cr y 

T] 6 



aH_y^ 

ri 4 ' 



aH 
~2^- 



(40) 



The problem for perturbations of the stream function 
and z-component of the vorticity vector reads 

VVi^-'^i, V2wi = 0, (41) 
which can be solved by applying boundary conditions 



dy 

dip 



i(x,0) = b{x) 



aH 



uJi{x, 0) 




(42) 



and an extra condition that reflects our definition of the 
stream function: 



^i(2;,0) = 
This can be solved to get 

Afo 2 



(43) 



ipi{x,y) = — —y + Pay 





y 


2 


An 




y 


2 


An 



e cos A„x 



f:(^P^ + ^^)e— sAnX, (44) 



n=l 

Conditions (Sll lead to 



y) = ^ + E (A^l'^e^"^ + M(2)e-^"^) cos(Anx). 

(45) 



p(l) - _ P(2) = _p 

Pn (-e 



T A 71 // _|_ p Xj^H 



^ 2Po 
H ' 

2Xr,He 



X„H\ 



and we obtain another dual scries problem, which is sim- 
ilar to (1291) and (1301): 



aoll + — j + ^an[l + 26A„F(A„H)]cos(A„a 



TT 

b—, 0<x<d/2, 
2ri 



(46) 



"0 + E cos(Ana;) = 0, (5/2 < .t < L/2, (47) 

n=l 

Here, 

cosh(2Ang) ~ 2\lH'^ - 1 ^ 
ao = Pq, CLn = ^ — — Pn. n > 1, 



and 



V{t) 



sinh(2i) - 2t 
cosh(2t) - 2^2 - 1 ■ 



(48) 



In the limit of a thin channel (where l^(t)|t_>oo 
2<"i + 0{t)), the dual series problem transforms to 



flo + 1 + 



3fe 



H + b 



• a„ cos(A„a;) = 

, 277 1 

n — 1 



2?7 1 + b/H' 
0<x< (5/2(49) 
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ao + ^ a„ cos(A„a-) = 0, 5/2<x<L/2, (50) 



which aUows one to evaluate 

aH bH(t)2 



an 



277 H + Ab- 3(j)2b 
The effective shp length is then j2ll | 

bHc^2 



H 



(51) 



(52) 



This exact equation represents a rigorous lower Wiener 
bound on the effective slip over all possible two-phase 
patterns in a thin channel. 

For completeness here we mention again the two lim- 
iting situations: 



4 01 



(53) 



(x) 

Interestingly, in this limit b^g does not depend on 6, be- 
ing a function of only H and a fraction of the gas area. At 

small b according to Eqs. (p4|) and ([54|). bl^ ~ b^g ~ b[g , 
so that the flow becomes isotropic. 

In the limit of a thick channel and sufficiently large 
local slip, we can simplify Eq.([2T|) and deflne the down- 
stream effective slip length as 



,11 



cos^ e + 6, 



cff 



(59) 



with b^'f^ given by Eqs. ([55]) and ([55)) . In the limit of 
perfect local slip it can be further simplified to get 



6iff(l+COs2e). 



(60) 



Using Eqs. ([55)) and ([5S)) we conclude that the flow is 

isotropic big ~ b^g 
the texture period. 



(x) 

b^g when b is much smaller than 



V. SIMULATION METHOD 



t>2 oc b. 



(54) 



In the thick channel limit, the dual series ([46)) and ([47)) 
take the same form as in prior work jsj (due to V{x — > 
oo) — > 1), whence we derive Q 



(55) 



The consideration as above of the same limits of small 
and large b give 



L 




In 


f n<f)2\ 






27r& 


sec(^— j+tan(^— j 



and 



"off l()<L<_f/ 



sec 



oc b 



TT(P2 

2 



(56) 



oc L, (57) 



C. Tilted stripes 



If the stripes are inclined at an angle 0, the effective 

(x) 

slip length of the channel, b^g can be calculated with 
Eq. ([2T)). provided that effective slip in eigendirections is 
determined from the numerical solution of ()29p , ([50)) and 
([46)) . ((47)). Some simple analytical results are possible in 
the limit of thin and thick channels. 

In case of a thin channel and large local slip, H <C 
min{&, L}, substitution of Eqs. ([33)) and ([53)) into Eq. 
([2T)) gives 



/'2 + 01 + 



301 cos^ 6 
ii — 302 cos-^ Q ' 



(58) 



To simulate fluid flow between parallel patterned plates 
a number of simulation methods could be used. These in- 
clude Molecular Dynamics, Dissipative Particle Dynam- 
ics, Stochastic Rotation Dynamics, classical Finite El- 
ement or Finite Volume solvers as well the LB method. 
Since the current paper does not address molecular inter- 
actions or liquid-gas transitions close to the surface, one 
can limit the required computational effort by applying 
a continuum solver for the Stokes equation together with 
appropriate boundary conditions to model local slip. As 
detailed further below, in particular in the thin channel 
limit with H <^ min{&, L} a very high resolution of the 
flow field is needed in order to measure 6off with required 
precision. Thus, simulation methods which require time 
averaging of the fiow field or finite lengths scales to re- 
solve a slip boundary render less efficient for the cur- 
rent problem. While Finite Element or Finite Volume 
solvers would be suitable alternatives, we apply the LB 
method 
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The LB approach is based on the Boltzmann kinetic 
equation 

'd_ 

di 



U • Vr 



(61) 



which describes the evolution of the single particle prob- 
ability density /(r, u,i), where r is the position, u the 
velocity, and t the time. The derivatives on the left- 
hand side represent propagation of particles in phase 
space whereas the collision operator $1 takes into account 
molecular collisions. 

In the LB method the time t, the position r, and the 
velocity u are discretized. In units of the lattice con- 
stant Ax and the time step At this leads to a discretized 
version of Eq. ()61l) : 



/fc(r + Cfc,t + l)-/fe(r,i) A: = 0,1, 



(62) 
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Our simulations are performed on a three dimensional 
lattice with B = 19 discrete velocities (the so-called 
D3Q19 model). With a proper choice of the discretized 
collision operator ft it can be shown that the flow behav- 
ior follows the Navicr-Stokcs equation [33. We choose 
the Bhatnagar-Gross-Krook (BGK) form [S^ 

f^A. = --(/fe(r,0-/r(u(r,t),p(r,<))) , (63) 



which assumes a relaxation towards a discretized local 
Maxwell-Boltzmann distribution /^"^ . Here, t is the mean 
collision time that determines the kinematic viscosity i> ~ 
'^^g^ of the fluid. In this study it is kept constant at 
T = 1.0. 

Physical properties of the simulated fluid are given by 
the stochastical moments of the distribution function. Of 
special interest are the conserved quantities, namely the 
fluid density p{r,t) = PoY^f. fk{'r,t) and the momentum 
p(r, t)u{r, t) = po J2k Cfc/fe(i", t), with po being a reference 
density. 

Within the LB method a common approach to describe 
the interaction between hydrophobic surfaces and the 
fluid is by means of a repulsive force term [3: [H; l39l - l4l| . 
This force applied at the boundary can be linked to the 
contact angle to quantitatively describe the wettability 
of materials [l^ |43 - |4^ . Alternatively, slip can be in- 
troduced by generalizing the no-slip bounce-back bound- 
ary conditions in order to allow specular reflections with 
a given probability (45l - l47| . or to apply diffuse scatter- 
ing [iSroOl l . The method we apply here follows the latter 
idea and uses a second order accurate on-site fixed veloc- 
ity boundary condition to simulate wall slippage. The on 
site velocity boundary condition is used to set a required 
slip length on the patterned surface. For the details of 
the implementation we refer the reader to [5l|, [s^ . 

Our geometry of configuration is the same as sketched 
in Figurc[l] but in simulations we employ periodic bound- 
ary conditions in x and z-direction, which allows to re- 
duce the simulation domain to a pseudo-2D system. To 
drive the fiow a constant pressure gradient is applied in 
X direction by means of a homogeneous acceleration, g, 
in the whole fluid domain. Even though the simulation 
domain can be reduced to be pseudo-2D, we find that the 
simulation of fluid flow with a large slip in the thin chan- 
nel limit still requires a system of several million cells 
in order to properly resolve the velocity field. The key 
issue is here that in order to resolve large slip lengths a 
minimum channel height is necessary. To reach the thin 
channel limit the stripe length has to be increased to even 
larger values. 

Following the definition given for the Fourier analysis 
in Sec. IIV[ all heights H and slip-lengths b are given 
non dimensionalised for a stripe length oi L = 27r. The 
resolution of the simulated system is then given by the 
lattice constant 



where J\f is the number of discretization points used to 
resolve the height of the channel. While systems with a 
height oi H = 1.0 can be simulated using a discretization 
of 1 X 32 X 200Ax^ only, decreasing iJ to 0.1 causes the 
required system size to be 1 x 70 x 4400 Aa;'^. To success- 
fully recover the exact results in the thin channel limit 
(H = 0.01) a system of size 1 x 104 x 64000Ax3, has to 
be simulated. 

The number of timesteps required to reach a steady 
state depends on the channel height, the velocity of the 
flow as determined by the driving acceleration as well as 
the fraction of slip and no slip area at the surface. For the 
simulations conducted in the thin channel limit a steady 
state velocity field exactly fitting the theoretical predic- 
tion develops after one to four million timesteps. In the 
thick channel limit, however, the number of timesteps re- 
quired can be an order of magnitude larger limiting the 
maximum feasible system height. Moreover, the transi- 
tion between slip and no slip stripes induces a distortion 
of the flow fleld with a range of ~ 3Ax. In order to keep 
the induced error below an acceptable limit, a minimum 
resolution of the channel length of 64Aa; is maintained. 
Additionally, the maximum flow velocity is limited due 
to the low Mach number assumption of our lattice Boltz- 
mann implementation. In effect, the acceleration mod- 
elling the pressure gradient has to be reduced increasing 
the time required for convergence. For example, a sim- 
ulation domain of 1 x 1024 x 64 Ax"^ as used to model 
the thick channel limit at i7 ~ 100 requires 12.5 million 
timesteps to equilibrate. 




Ax = 



HL 
2^' 



(64) 



FIG. 2: Illustration of the measurement of effective slip eigen- 
values from lattice Boltzmann simulations. The velocity infor- 
mation of the whole domain is projected onto a single plane. 
To this cloud of data, Eq. (|65p is fitted, effectively averaging 
the flow fleld over the whole channel. 



We compare measurements of 6eff by permeability esti- 
mates and velocity profiles, respectively. The permeabil- 
ity is calculated from measurements of the flow rate, ac- 
cording to Eq. . This allows to determine the effective 
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slip by Eq. p^ . Alternatively, the profiles of the velocity 
in flow direction are averaged over the whole system by 
projecting the velocity information of the whole domain 
onto a single plane. Then, the effective slip length bes 
is found by a Levcnbcrg-Marquardt fit of the adjusted 
Hagen-Poiscuillc equation 



Uxiy) 



gH^ ({y-Hf {y - H) 



2-q 



{H + 6eff ) 



(65) 



with the known g and H (see Fig. [5] for an illustration). 

The error of the effective slip measurements is de- 
termined by two factors, namely the resolution of the 
channel height and the absolute slip length of the par- 
tially slipping stripes. For poorly resolved channels with 
10 < A/" < 30 and small slip lengths the permeability 
measurements still produce accurate results, whereas a 
fit of the velocity profiles fails. If 6 ^ iJ, however, the 
quality of the obtained data declines. For an increase 
in resolution (30 < J\f < 100) both approaches allow 
very precise measurements for intermediate slip lengths 
of up to two orders of magnitude larger than the chan- 
nel height. However, if b is increased further, due to 
discretization the error in the effective slip determined 
by the permeability measurement increases significantly 
rendering this method inefficient since a higher resolu- 
tion would be required. For example, to keep the error 
in the determination of a prescribed slip in the order of 
10^ Aa; below 5% the permeability method requires the 
channel height to be resolved by 200 lattice sites, while 
for the measurement by fitting the velocity profiles 100 
sites suffice. 

To validate the concept of a tensorial slip by simula- 
tions of a flow past tilted stripes, we do not rotate the 
surface pattern with respect to the lattice, but instead 
change the direction of the acceleration in the ?/ — z-plane. 
This avoids discretization errors due to the underlying 
regular lattice occurring in case of a rotated surface pat- 
tern. We extract the downstream slip by projecting the 
slip measured on the main axes onto the pressure gradi- 
ent direction. 



VI. RESULTS AND DISCUSSION 

In this section we compare results of our LB simula- 
tion with analytical example calculations and numerical 
solutions of the dual series 1^ and (gB]), (gT]) (sec 
Appendix . 

As a benchmark for the simulation, we start with a 
thin channel, where striped surfaces were shown to pro- 
vide rigorous upper and lower Wiener bounds on the ef- 
fective slip over all possible two-phase patterns [2l|. In 
order to reach the thin channel limit, a dimcnsionlcss 
height oi H ^ 0.01 is chosen. The slip lengths are set to 
6 = 10-^ H {0.1 Ax) and b = 10^ H (102000Aa;), differing 
each three orders of magnitude from the channel height 
and reaching the limits of small- (cf. Eqs. and ([M)) ) 
as well as large slip (Eqs. and ([55)) ). As preparatory 




FIG. 3: Eigenvalues of the effective slip length tensors sim- 
ulated in the limit of a thin channel (symbols). The lines 
represent results of theoretical calculations by Eqs. p2|l and 
(|52|) . The data show that at small b the eigenvalues of bcff 
decrease as compared to a large local slip at the gas sector, 
and that the slip-length tensor becomes isotropic resulting in 
b^g to become hardly distinguishable in the 6 = 10 H case. 



tests showed, a minimum channel height oi H = lOOAx is 
required to measure slip-lengths of 6 = 10^ H correspond- 
ing to 10^ Ax. For a dimensionless height oi H ~ 0.01, 
we choose a simulation domain of 1 x 104 x 64000Aa;'^. 
For each of the two slip-lengths, longitudinal and trans- 
verse flow was simulated for a different fraction of sur- 
face gas phase, ranging from no-slip (02 = 0) to homo- 
geneous partial slip ((/)2 = 1). The local acceleration 
here (and below for thin channel simulations) was kept 
at g = IQ-^Ax/At^. 

Fig. [3] shows the exact eigenvalues of the effective slip 
tensor in the thin channel limit, Eqs. ([32]). ([52)) . for both 
slip lengths b. The flt of the simulation data and the 
analytical limits is excellent for all separations. In the 
case of small local slip in the thin channel the effective 
slip remains isotropic despite of the inhomogeneity of the 
boundary. For large local slip, we observe truly tensorial 
effective slip and highly anisotropic flow over the surface. 
These simulations demonstrate that flnite size effects and 
resolution effects are well controlled, and the size of the 
system is sufficient to avoid artifacts. Another impor- 
tant point to note is that in our theoretical analysis all 
equations were derived ignoring stripe edge effects. An 
excellent agreement between theoretical and simulation 
results indicates that the edge effects do not influence the 
simulation results signiflcantly. 

Fig. ID shows the eigenvalues of the effective slip length 
tensor as a function of 02 for a thick gap. For these 
simulations the acceleration has been reduced down g ~ 
10~^Aa;/At^ to obey the low Mach number (see Section 
W\ limit. The time to reach a stable state increased then 
to 15 • 10^ At. Simulation results are presented for two 
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FIG. 4: Simulation results of eigenvalues of bcfi as a function 
of fraction of gas sectors, (j>2, in the limit of thick channel 
(symbols). Lines represent corresponding theoretical values 
obtained by a numerical solution of ((29)) . pO)) and (|46)) . (|47)) . 



different slip lengths of 5 = 1.0 and b = 10.0 in a sys- 
tem oi H = 0.1, where L is now resolved by 4400 lattice 
sites. The theoretical solutions represented by the lines 
were obtained by the dual series approach. We find that 
the fit is excellent for all fractions of the slipping area, 
indicating that our semi-analytical theory is extremely 
accurate. The data presented in Fig. U show larger ef- 
fective slip for a lower slip to height ratio, i.e. a thicker 
channel. This illustrates well the earlier suggestion that 
effective boundary conditions for this channel geometry 
are controlled by the smallest length scale of the prob- 
lem m. 

To check the validity of the tensorial slip approach, we 
now orient the texture relative to the x-axis, which in our 
model is always aligned with the applied pressure gradi- 
ent. Fig. [S] and [5] show two sets of effective downstream 
slip lengths simulated with several 8, but fixed H — 0.1 
and (/)2 = 0.5, which results in a maximum transverse 
flow in a thin channel situation [22| . 

In the first set (Fig. [5]), we consider thin channels and 
vary b/H from 1 to 1000. Fig. [5ji shows simulation data 
obtained using a channel of height H = 0.1. Further, 
theoretical curves calculated with Eq. [T] arc presented. 
Here, eigenvalues of the slip-length tensor are obtained 
by numerical solution of the dual series. The fits of the 
simulation data are in very good agreement with the nu- 
merical solutions of the dual series suggesting the validity 
of the concept of a tensorial slip in a thin channel situa- 
tion. Note that the simulation results of Fig. [5^ cannot 
be compared with the analytical expression, Eq. (j58p , be- 
cause Fig. \^ is based on a relatively moderate value of 
local slip at the gas sectors, whereas Eq. ([55|) requires very 
large b. To validate predictions of this analytical formula, 
the channel height was decreased down to H = 0.01. 
Simulation results are presented in Fig. [SJd, confirming 




FIG. 5: Downstream effective slip lengths of the thin channel 
simulated at <j}2 — 0.5 for stripes inclined at different angle 
O (symbols). Especially in the limiting case (b) the value of 
has to be calculated from the measured eigenvalues by 
Eq. (|21|I . All the lines are predicted theoretically downstream 
slip lengths: (a) Calculated by using Eq.([T]) [or Eq.((2T])] with 
eigenvalues of bcfi determined from numerical solutions of 
(ESI),!® and (ge]),!!?]); (b) Calculated with Eq.([58l). 



the surprising accuracy of a simple analytical expression, 
Eq. (|58p . We remark that in this important limit of valid- 
ity of Eq. ((58|) . b^^ /H is quite large, although 6^^^ itself 
is small. This may have implications for a reduction of 
a hydrodynamic drag force |23l . [23 | . Also passive mixing 
might be an interesting application, since the anisotropy 
of flow is very large, which is optimal for a transverse flow 
generation [22|. We suggest that our simple asymptotic 
result could be intensively used to simplify theoretical 
analysis of these important phenomena. 

In the second set as shown in Fig. [H] a thick channel 
(of height H = 100) is simulated. Fig. [S^ plots simula- 
tion results for several b/H varying from 10-3 to 1( sym- 
bols). Similarly to previous examples, we found a very 
good agreement between simulation results and predic- 
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FIG. 7: Effective downstream slip lengths at different channel 
thickness {<f)2 = 0.75, b/L = 5.0). Symbols illustrate the 
simulation data, and curves show theoretical predictions. 



ary condition is not a characteristic of the liquid-solid in- 
terface alone, but depends on the flow configuration and 
interplay between the typical length scales of the prob- 
lem. Again, the simulation and theoretical data are in 
the excellent agreement for longitudinal, transverse, and 
tilted stripes. Thus, Fig. [7] unambiguously shows that 
the tensorial slip boundary condition, originally justified 
for a thick channel, to any channel thickness can be gen- 
eralized to any channel thickness. 



FIG. 6: Effective downstream slip lengths for tilted stripes 
(the thick channel limit) simulated at (/!>2 = 0.5 (symbols). All 
the lines are predicted theoretically downstream slip lengths: 
(a) Calculated by using Eq.(IT]) [or Eq.((2T|)] with eigenvalues 
of bcfi determined from numerical solutions of (|29|l , (|30p and 
(SSJ, (SZl); (b) Calculated with Eq. ((59)) with eigenvalues eval- 
uated with Eqs. ^7} and ((57)) . 



tions of tensorial Eq. ([T]) with eigenvalues computed with 
our semi-analytical theory. We found that already for 
the case b = H ^ 100 our simulations reach the limit 
of large slip in the thick channel, so that the compar- 
ison with analytical solutions is possible. To examine 
this more closely, the simulations results obtained in this 
limit are reproduced in Fig. |B]d. Also included are the 
theoretical results calculated with asymptotic formulae, 
Eqs. ([57]) and ([57]) . which perfectly fit simulation data. 

Finally, we simulate the downstream slip length as a 
function of the channel thickness with the focus on the 
intermediate gap situation. Fig. [7] shows the typical sim- 
ulation results (the example corresponds to b/L = 5.0 
and 02 = 0.75) and demonstrates that the effective slip 
lengths increase with H and saturate for a thick gap. 
This fully confirms the statement that an effective bound- 



VII. CONCLUSION 

We have investigated pressure-driven flow in a flat- 
parallel channel with one hydrophilic and one super- 
hydrophobic surface, and have given some general the- 
oretical arguments showing that a concept of an effec- 
tive tensorial slip is valid for any thickness (compared 
to a super-hydrophobic texture scale). The eigenvalues 
of the effective slip-length tensor depend on the gap, so 
that they cannot be viewed as a local property of the 
super-hydrophobic surface, except in the thick channel 
limit. Instead, the slip-length tensor represents a global 
characteristic of the channel. The mathematical proper- 
ties of the slip-length and permeability tensors allowed 
us to derive a simple analytical formula for an effective 
downstream slip length in case of inclined to a pressure 
gradient textures. Our analysis is validated by means of 
LB simulations. 
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cos(mz) and Eq. IA6I by cos(mi) and we then integrate 
over [c, tt] . The resulting equations are summarized to 
obtain a system of linear algebraic equations, 



(A7) 



which can be solved with respect to a„ by standard nu- 
merical algebra tools. 

The coefRcients for the longitudinal case are (m > 0) 



Appendix A: Numerical method 



Equations ^ and (|46j), gZ]) provide a complete 
description of hydrodynamic flow and effective slip in 
eigendirections. Their exact solution is possible for the 
limits of a thin and a thick channel only. In order to solve 
the problem for general channel thickness, the following 
numerical algorithm has been used. 

It is convenient to change to dimensionless values. We, 
therefore, choose L / {2tt) as the reference length scale and 
aHL I {4:Trr]) as the velocity scale. We make the substitu- 
tion 

{x,y,z) ^ -^{x,y,z), H^-^H, b^-^b, (Al) 
Ztt ztt Ztt 



aHL. 



n > 0, 



(A2) 



where non-dimensional variables are denoted by tildes. 
This procedure gives the dual series problem for longitu- 
dinal flow in the form 




zsm{mz)dz+ / cos{mz)dz, (A8) 



,11 l + bncoth{nH) f . ^ . , 

Al^m = / sm(TOz) smlnzjdz 

n J 





-|- / cos(mz) cos(7iz)(i2, n>\, 



(A9) 



Bl^ = b J 5sin(mz)dz. (AlO) 





For transverse flow we get 

_L _ 1 + 2bnV{nH) 



sin(mz) s'm{nz)dz 



.0 1 + ^ 



dn 1^1 + 6ncoth(ri_ff) 

n=l 

&, < z < c, 



(A3) 



do + ^ dn cos(nz) =0, c < z < tt, (A4) 

n=l 

where c = tt4>2 = tt6/L. Similarly, the equations for the 
flow in the direction orthogonal to the stripes is written 
as 



do 1^1 + j + dn [l + 2bnV{nH) cos( 
= 6, < i; < c. 



nx 



(A5) 



-I- / cos(mz) cos(nz)dz, n > 1 (All) 



(A12) 



4^ — 4" — rII 



For the numerical evaluation, the linear system is trun- 
cated and reduced to a. N x N matrix and the solution 
is then found to converge upon truncation refinement. 
According to the definition 



(Us) 



(A13) 



the dimensionless and dimensional effective slip lengths 
are given by 



do + 2. o-n cos{nx) =0, c < i < tt. (A6) 



n=l 



After integrating Eq. IA3I over [0, z] and Eq. IA5I over 
[0,i]), we multiply the result by sin(mS) (sin(TOi), re- 
spectively), where m is a nonnegative integer. We then 
integrate again over [0,c]. Eq. IA4I is multiplied by 



^off 



1 - do/H 



Ztt 



(A14) 



(A15) 
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